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ABSTRACT 

We report on the observation of a very small glitch observed for the first time in a millisecond pulsar, 
PSR B1821-24 located in the globular cluster M28. Timing observations were mainly conducted with 
the Nangay radiotelescope (France) and confirmation comes from the 140ft radiotelescope at Green 
Bank and the new Green Bank Telescope data. This event is characterized by a rotation frequency step 
of 3 nHz, or 10 -11 in fractional frequency change along with a short duration limited to a few days or 
a week. A marginally significant frequency derivative step was also found. This glitch follows the main 
characteristics of those in the slow period pulsars, but is two orders of magnitude smaller than the 
smallest ever recorded. Such an event must be very rare for millisecond pulsars since no other glitches 
have been detected when the cumulated number of years of millisecond pulsar timing observations 
up to 2001 is around 500 for all these objects. However, pulsar PSR B1821-24 is one of the youngest 
among the old recycled ones and there is likely a correlation between age, or a related parameter, and 
timing noise. While this event happens on a much smaller scale, the required adjustment of the star 
to a new equilibrium figure as it spins down is a likely common cause for all glitches. 
Subject headings: stars: neutron - stars: rotation - pulsars: general - pulsars: individual(PSR B1821- 
24) 



1. INTRODUCTION 

In long period pulsars rotational irregularities are ob- 
served that are thought to result from erratic internal 
structure changes. Two types of irregularities are ob- 
served: "timing noise", which is a slow unpredictable 
fluctuation of rotation; and "glitches" , which are sudden 
changes in the rotation rate. Glitches mainly occur in 
younger pulsars (Lyne, Shemar & Graham-Smith 2000), 
and are usually interpreted as a sudden transfer of angu- 
lar momentum from a faster-rotating component of su- 
perfluid to the solid crust of the star (Baym et al. 1969, 
Anderson & Itoh 1975). These phenomena provide an 
opportunity to study the internal structure of neutron 
stars. Days or months are necessary to recover after the 
glitch and such long relaxation have been interpreted as 
strong evidence for superfluid in the neutron star (Sauls 
1989). Glitches have been observed with fractional fre- 
quency step amplitudes down to the level of 10~ 9 by She- 
mar & Lyne (1996). They further estimated that 90% 
of the glitches with size greater than 5 x 10 -9 are de- 
tected. Despite carefull inspection of many years of high 
precision timing results of millisecond pulsars by several 
groups, no glitches have been observed in this subcate- 
gory of pulsars. This fact is consistent with the statis- 
tical result that the oldest pulsars, like the millisecond 
pulsars, are relatively free from glitches (Lyne, Shemar & 
Graham-Smith 2000). The Nangay radiotelescope is in- 
volved in a long term timing program on two millisecond 
pulsars for more than 15 years. In this Letter, we report 
the first micro-glitch observed in a recycled millisecond 
pulsar which is two orders of magnitude lower than any 
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previously detected. 

2. OBSERVATIONS 

The timing observations presented here have been con- 
ducted at the large decimetric radio telescope (7000 
square-meters, equivalent to a 93m dish) located near 
Nangay (France). We have gathered observations of 
PSR1821-24 since 1989 October about 10 times per 
month at 1.4 GHz. At Nangay, the pulsar signal is co- 
herently de-dispersed prior to detection by using a swept 
frequency local oscillator (first at 80 MHz and then at 150 
MHz since 2000) in the IF chain. The pulse spectra are 
produced by the station digital autocorrelator. While 
most timing observations use a fixed center frequency 
and determine a pulse arrival time at this frequency, our 
swept-LO approach leads us to determine a frequency of 
arrival for the start time of the observation by a stan- 
dard cross-correlation of the daily integrated pulse with 
a pulse template (Cognard et al. 1996). The Nangay ob- 
servations, which are occasionally done at 1.68 GHz to 
monitor the variations of the column density of dispers- 
ing electrons, were suspended between 1999 November 
and 2000 July owing to the upgrade of the telescope. 
Observations of B1821-24 were also conducted with the 
140ft Green Bank telescope between 1989 and 1997, and 
briefly with the new 100m Green Bank Telescope (GBT) 
during 2002 and 2003. 

The analysis of the pulsar timing data has been car- 
ried out by our software package (AnTiOPE) that fits the 
observations to selected parameters involving the pulsar 
spin and astrometry along with delays due to disper- 
sive effects (Cognard et al. 1995). The Earth's orbital 
motion is obtained from the Jet Propulsion Laboratory 
ephemerides DE407 and the adopted timescale is Inter- 
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Fig. 1. — Timing residuals of PSR1821-24 at Nangay Radio Observatory from October 1989 to August 2003 around radio frequency 
1.4 GHz. Crux markers from 1989 to 1997 and in late 2002 and early 2003 are Green Bank TOAs. A set of parameters (period, period 
derivative, position, proper motion plus few minors ones) is adjusted over the Oct 1989 to Feb 2001 period and used to minimize the TOA 
residuals. After March 2001, both TOAs from original set of parameters and TOAs obtained with a <5P/P change of 10 -11 arc shown. 
The inset is a enlargement of the transition region showing that this transition does not last more than a few days or a week. Also clearly 
visible is the effect of the extra dispersive delay due to the occultation by the solar corona occuring each year around Christmas (Cognard 
et al., 1996). 



national Atomic Time (TAI). If a model for the Times 
Of Arrival (TOAs) is complete, then the residuals, ob- 
served TOAs minus model TOAs, will have a Gaussian 
distribution about zero. A few months after 2001 March, 
we found that the standard model of B1821-24 was in- 
adequate. Figure 1 shows the timing residuals based on 
a model up to 2001 March and the abrupt change in 
slope - a glitch in the rotation frequency. Various checks 
were done regarding the instrumentation. The most con- 
vincing evidence for a real effect came from the lack of 
any similar effect in the analysis of timing observations 
of other millisecond pulsars (MSPs) in the Nangay pro- 
gram. We also checked that the effect was not a function 
of the observing radio frequency and therefore could not 
be the result of dispersion measure variations (Cognard 
& Lestrade 1997). Moreover, the Green Bank measure- 
ments, with two different telescopes before and after 2001 
March but at the same frequency, provide confirming ev- 
idence (Fig. 1). 



3. RESULTS 



The parameters of a glitch in rotation (phase step, per- 
manent frequency step, permanent frequency derivative 
step) were added to AnTiOPE so that we could obtain 
a consistent overall fit. Table 1 shows the different fitted 
parameters. We tried to fit for a phase step at the date of 
the glitch, but nothing reliable was found. A step of the 
order of 15 (is was found, but at the price of a very clear 
discontinuity in the TOA residuals. We decided not to 
include this parameter in the fit. A frequency derivative 
step was also found, but its significance was marginal 
and the presence of timing noise is a significant source of 
systematic error. The short duration of the event (inset 
of Fig.l), which is limited to less than a week, is the most 
important fact that leads to our glitch interpretation, as 
opposed to timing noise as discussed below. 

Because of the dramatic change seen, we are driven to 
also explore the pulsar spin frequency variability by bin- 
ning data before and after the glitch. A robust approach 
was then chosen with independent fits over different sub- 
sets of data. Fits for frequency (and arbitrary rotational 
phase) were done over data sets with durations of 100 
and 365 days. The position and proper motion param- 
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Fig. 2.— Evolution of the PSR B1821-24 rotational frequency. 
Variations around a mean value are plotted in nanoHertz as a func- 
tion of the date. Black points and line are for frequencies deter- 
mined over a full year (365days) while gray points were obtained 
over shorter data spans of only lOOdays. Horizontal error bars in- 
dicate the duration of the box used to derive the value. The glitch 
occuring in March 2001 is clearly visible as a step of the order of 
3nHz. 



eters from the whole data set were used to fix the posi- 
tion in each data subset. The results in Figure 2 show a 
frequency step of approximately 3 nHz, or 10 _ 11 in frac- 
tional frequency change, occuring in March 2001. Note 
that in both cases, 100 days or 365 days, the data bins 
were phased with the date of 2001 March 10 when the 
event is assumed to have occured. Despite some obvi- 
ous noise in the rotational frequency (e.g., during 1996) 
the frequency step is confirmed. The frequency residu- 
als could be further characterized as showing a glitch in 
frequency derivative in late 1996 with magnitude of 0.5 
nHz y _1 , or 1.7 x 10~ 17 s -2 ; we will not discuss this 
possibility further in this short account. However, note 
the long term nature of this episode contrasting with the 
sudden frequency change of the 2001 glitch. 

The 2001 March frequency glitch is two orders of mag- 
nitudes lower than any previously observed. The small- 
est ones were observed at the level of 10 -9 in pulsars: 
0525+21 (Downs 1982), 1907+00 (Gullahorn ct al. 1976) 
and 0740-28 (D'Alessandro et al. 1993). However dis- 
crimination between glitches and timing noise is difficult 
at this level. Shemar & Lyne (1996) estimated that 90% 
of the glitches with relative frequency step greater than 
5 x 10 -9 were detected. Here, we are observing a millisec- 
ond pulsar with a typical TOA uncertainty of 2 /us, which 
together with a high timing precision and stability allows 
easy detection of a glitch as small as 10~ n . This small 
event may indicate a new regime in the evolution of neu- 
tron star internal structure. Indeed there are some very 
small glitches which may have been identified of both 
signs in term of rotational frequency increment from the 
very young Vela pulsar (Cordes et al. 1988). More events 
in B1821-24 and other MSPs are required to provide a 
clearer picture. Part of the surprise and signifiance is 
that other MSPs are so stable. PSR B1937+21 is known 
to be unstable at the microseconds level (Kaspi et al. 
1994, Lommen & Backer 2001), but no glitches down to 
the 10 -11 level have been seen in over 20 years. The 
same is true for many other MSPs. We estimate that the 



TABLE 1 

PSR B1821-24 PARAMETERS OBTAINED WITH ANALYSIS SOFTWARE 
ANTlOPE, TAI TIME SCALE AND JPL-DE407 EPHEMERIS. 



Spin and astromctry parameters 

Frequency (Hz) 327.405641011622(5) 

Frequency 1st derivative (Hz/s) -173.5288(l)e-15 
Frequency second derivative (Hz/s2) -4.15(7)e-26 

Epoch 2451468.5 

Right Ascension J2000 18 24 32.00821(3) 

Declination J2000 -24 52 10.784(8) 

RA proper motion (mas/yr) 0.27(8) 

DEC proper motion (mas/yr) -15.2(16) 
Glitch parameters 

Glitch Epoch 2451980.0 

Glitch Frequency step (Hz) 3.12(3)e-9 
Glitch Frequency derivative step (Hz/s) 5.9(5)e-18 



cumulated number of years of millisecond pulsar timing 
observations up to 2001 is around 500 for all objects. 

In order to assess the amount of potential timing noise 
present in B 182 1-24 data, we performed a fit of harmoni- 
cally related sinusoids to the timing residuals in the inter- 
val preceding the 2001 March micro-glitch following the 
work of Hobbs (2002). The rms of the curve produced 
can be seen as a measure of the timing noise present in 
the data. The value obtained for B1821-24 is 24 fj,s. This 
value is a reasonable extrapolation of the trend of rms 
vs period derivative P found by Hobbs (2002) and ear- 
lier authors. The value is larger than but comparable 
to that of the 1.5ms pulsar B1937+21, but significantly 
larger than that for other MSPs (Backer 2004). 

In the known large glitches in slow pulsars, the sud- 
den increase in rotational frequency is followed by an ex- 
ponential recovery and a possible permanent change in 
period derivative. Glitches with exponential recovery ob- 
served in young pulsars indicate that the neutron star is 
probably made of two components, crust and supcrfluid 
interior, that are partially coupled together. A fraction 
of the quantized, superfluid vortices can be pinned to 
defects in the crust to trap angular momentum as the 
star spins down and the vortices migrate outward (Alpar 
et al. 1981). The trapped angular momentum is sud- 
denly released at the time of a glitch by an unexplained 
catastrophic unpinning event where angular momentum 
is transfered from the inner part of the crust to the crust 
to produce the observed speed-up. 

Glitch activity is thought to be maximum for young or 
adolescent pulsars, with ages up to 2 x 10 4 y. Ages of 
pulsars are usually estimated from the rotation evolution 
with a characteristic age given by P/2P assuming a mag- 
netic dipole braking index of 3 (Lyne & Graham-Smith 
1998). Determination of the age of a pulsar in a globular 
cluster can be corrupted by the gravitational potential 
through the period derivatives measurements. However, 
Phinney (1993) showed that contribution from the clus- 
ter to the period first derivative is limited to 6% with 
a P term of 0.3 -16 s _1 at maximum. The characteris- 
tic age of 30 x 10 6 years for B 182 1-24 is therefore fairly 
well established, and makes this millisecond pulsar the 
youngest among the globular cluster and galactic field 
MSPs. 

Cordes et al. (1988) find a bimodal distribution of 



discontinuities in the timing record of the Vela pulsar. 
The large events are the conventional "glitches" showing 
sudden spinup and subsequent exponential decay, and 
amplitudes (|Az//z/|) of ~ 1CP 6 . They call the small 
events "micro jumps" . These have amplitudes of ~ 10~ 9 
and are not merely scaled down versions of the glitches. 
They suggest that random timing noise in other pulsars 
is caused by unresolved microjumps. With just a sin- 
gle event in the ensemble of millisecond pulsar timing 
records, it is difficult to speculate on how the glitch, jump 
and timing noise properties evolve with age and mag- 
netic field. For example, Lyne et al. (2000) say "there is 
a broad range of glitch amplitude, and there is a strong 
indication that pulsars with large magnetic fields suffer 
many small glitches while others show a smaller num- 
ber of large glitches." Following this to millisecond pul- 
sars with small magnetic fields, they should exhibit small 
number of large glitches. One microglitch does not fit in 
with this statement beyond suggesting that a wide distri- 
bution of glitch amplitudes is possible. Glitch amplitude 
distribution and timing noise (e.g., Backer 2004, Hobbs 
2002) in the millisecond pulsar population will give us 
more insight into the nature of neutron star interiors in 
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the future. 



4. SUMMARY 

For the first time, a glitch was observed in a recycled 
millisecond pulsar. We must note that this pulsar ap- 
pears to be one of the younger MSPs. Long term ob- 
servations of recycled pulsars are needed to assess the 
effect of this behavior and influence on long term studies 
like the search for background of primordial gravitational 
waves. If other glitches are detected in very stable recy- 
cled pulsars, the perspective of such long term goals will 
need to be thought again. 
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